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Abstract 

Density functional theory is discussed in the context of one-particle systems. We show that the 
ground state density po(x) and energy Eq are simply related to a family of external potential energy 
functions with ground state wave functions t^ n (x) oc po(x) n and energies E n = 2nE$ for certain 
integer values of n. 
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I. INTRODUCTION 



The objective of density functional theoryA (DFT) is to transform the many-particle prob- 
lem into an equivalent one-particle problem. Therefore, it seems pedagogically reasonable 
to give an introductory level presentation of DFT, in the context of one-particle systems, 
that (briefly) ovoids references to the complications imposed by systems with two or more 
particles. Several papers^ devoted to DFT have appeared in this journal. 

The basic idea of DFT is to transform from the ground state wave function ipn to the 
one-particle density function p as a focus of attention. This has obvious advantages for 
many-particle systems. For one-particle systems, it allows the exact ground state energy E 
to be determined via the variational principle with respect to the one-particle density. For 
one-particle systems 

Po(x) = \^ (x)\ 2 , (1) 
where the density function satisfies the normalization condition 



p (x)dx = 1. 

The wave function ^ satisfies the time independent Schrodinger equation 

h 2 d 2 



+ V(x) 



(2) 



(3) 



2m dx 2 

where V(x) is the external potential energy function for the particle of mass m. We consider 
systems where ipo is real, so that the energy may be written as 

h 2 d 2 



En 



Po{x) 



2m dx 2 

We now introduce the energy density functional 



+ V(x) 



po(x)dx. 



E[p) = T[p] + J p(x)V(x)dx, 



where the kinetic energy density functional T[p] is 

h 2 d 2 



T[p) 



p{x) 



p(x)dx. 



(4) 



(5) 



(6) 



2m dx 2 

The arbitrary density function p satisfies the same normalization condition as pn in Eq. 0. 
According to the Hohenberg-Kohn 1 (HK) theorem, (1) the external potential V(x) is deter- 
mined uniquely 3 by po(x), and (2) the energy density functional E[p] satisfies the condition 



E[p] > E 



(7) 



Variationally speaking, the minimum value of E[p] with respect to p yields the ground state 
energy Eq = E[po}. The density functional T[p] may be written more conveniently as 

1 dp(x) 



p(x)dx (8) 



p{x) dx 

According to the variational principle 

5 (e[ P ) -xj p(x)dx^j = 0, (9) 

where A is a Lagrange multiplier. In terms of variational (or functional) derivatives^ 

8E[p] 



8p(x) 

so that the equation for the density becomes 



A, (10) 



£(^ ln "W 2+ (s ln ' , w)V^) = A. (11) 



This equation may be rewritten in a form similar to Eq. (J3J): 

h 2 d 2 



2m dx 2 ' 



p(x) = 2Xp(x), (12) 



where the effective potential V e g is 



v^M)=m X ) + ^(~P^) . d3) 

Equation ()12|1 is an example of the nonlinear Schrodinger equation that is encountered in 
DFT calculations. The standard method of solution is to solve it iteratively, starting with 
an initial guess for V e g- At the solution point p = po and A = Eq. Equation (fTTj) may be 
transformed into a first order differential equation for y{x) = dlnp(x)/dx: 

-H 2 t^ + ^} +v(x)=x - (i4) 

This is a version of the Riccati 5 equation. A solution method for this equation has been 
discussed in some detail by Haley. 6 

In Sec. [Hi we derive equations that give the external potentials and ground state energies 
for the family of ground state wave functions ip n {x) oc po(x) n , where n = 2 J and j (j > 0) 
is an integer. In Sec. IIII[ we apply the results of Sec. |H] to several well known systems. In 
Sec. IIV1 we give a summary and some closing remarks. 



II. GENERATING EXTERNAL POTENTIALS WITH THE GROUND STATE 
DENSITY 



In this section we demonstrate one consequence of the HK theorem: The ground state 
density po(x) determines a family of external potentials and their respective ground state 
energies. We give some examples in the next section. 

Consider the wave functions given by 



i>n{x) = c n p {x) n , 



(15) 



where the integer n = 2 J for the integer j, and the constants c n are chosen to normalize 
ip n (x). Each ground state wave function given by Eq. (JT5|l satisfies the Schrodinger equation 

h 2 d 2 



2m dx'" 



+ VJx) 



ip n (x) = E n ijj n (y), 



(16) 



where E n is the ground state energy for the external potential V n (x). For example, applying 
Eqs. (|T2|) and (pH?|) iteratively for n — 1, 2, and 4 yields 



Vi(x) = V eS (x; [po]) 

h 2 ( l d 

= 2V(x) 



V 2 (x) = 2V 1 (x) + (2) 



4m \ po (x) dx 

h 2 ( 1 d 



4m \ po ( x ) dx 



V 4 (x) = 2V 2 (x) + (4) z — 



h 2 ( 1 d 



-po(x) 



(17) 
(18) 
(19) 



Am \po(x) dx 

In the Appendix, we show that the family of external potential energy functions and ground 
state energies are given, respectively, by 



V n (x) = 2nV(x) + (2n -l)n 



E„ = 2nE n 



h 2 ( 1 d 



Am \ p (x) dx 



po{x) 



(20) 
(21) 



III. EXAMPLES 

We have applied Eq. ([20)1 to (i) the particle in a box, (ii) the harmonic oscillator, and 
(Hi) the attractive delta potential. The results for these systems are summarized in Table UJ 
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IV. SUMMARY AND REMARKS 



We have presented DFT for one-particle systems. As a consequence of the HK theorem, 
the ground state density po(x) determines the external potential V(x) and a family of related 
external potentials V n (x) given by Eq. (|2(jp . Moreover, the ground state energies for these 
external potentials are given by Eq. (|21|) . These results results are directly applicable to 
higher dimensional system by replacing x by the position vector r and the derivative operator 
d/dx by the vector operator V. 

We hope that our presentation will be useful and insightful to teachers and students of 
quantum mechanics. 
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APPENDIX A: A DERIVATION OF EQUATIONS d20j) AND (HH) 
1. Equation (HSu]) 

The equations for the external potentials discussed in Sec. |H] may be transformed into 
the linear difference equation 

m^ 2 



z j+1 = 2z j +[V +L ) fix), (Al) 
where Zj = V 2 j(x) and the function f(x) is 

The solution 7 to Eq. (jAl|) is 

Zj = 2^ + 2V(x)E2 m 

= 2^ o + 2^70r)(2 J -l). (A3) 

In terms of n = 2 J 

V n (x) = nV 1 {x) +2n{n- l)f{x). (A4) 
Inserting the expression for V\{x) given by Eq. (fTT|) into this equation yields Eq. (pjj) . 
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2. Equation (f2"T|) 



Consider the two equations resulting from the insertion of the ground state densities po(x) 
and p n (x) oc p (x) 2n , respectively, into Eq. (JTTJ): 

~^{ 2 ' ! ^ lnMl)2+4 " 2 (^ lnpo(l) ) J + l «W = B »- ( A6 > 

Multiplying the first equation by 2n and then subtracting the resulting equation from the 
second equation yields Eq. (J5TJ). 
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TABLES 



TABLE I: The external potential V(x), the ground state density po(x), the ground state energy Eo, 
and the external potential V n {x) given by Eq. (120(1 for (i) the particle in a box, (H) the harmonic 
oscillator, and (Hi) the attractive delta potential. 



System 


V(x) 


Po(x) 


Eq 


V n (x) 


(0 


{ for \x\ < L/2 
j oo otherwise 


(2/L) cos 2 (ttx/L) 


(Hit/L) 2 /2m 


2n (2n - 1) E tan 2 (irx/L) 


(H) 


mu 2 x 2 /2 


yjmuj 1 'nh exp (—mujx 2 /h) 


hto/2 


4n 2 V(x) 


(Hi) 


-gS(x) 


{ing/h 2 ^ exp {—Img \x\ /h 2 ^j 


-mg 2 /2h 2 


2nV(x) - 2n(2n- 1) E 
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